This assignment is for ETC5521 Assignment 1 by Team goanna comprising of Chengzhi Ye, Emily Sheehan, and Dea Avega Editya.

1 Introduction and motivation

The Australian climate is generally hot and dry, which means that most regions can be affected by bushfires at anytime of the year (Australia 2020). Bushfires vary in their magnitude and temperature. Some bushfires can go on for days, weeks or even months. Some bushfires are out of control, while others can be contained.

Last year, the bushfires in Victoria and New South Wales captured the attention of people worldwide. They caused destruction and devastation for several months. Around 33 lives were lost, over 1 billion mammals died, and over 3,000 homes were destroyed (Lisa Richards, Nigel Brew, 2020).

This analysis hopes to understand the relationship between climactic conditions and bushfires, and determine whether climate change has influenced the number of bushfires.

R has been used as the main tool for cleaning and analysis. The analysis proceeds as follows; the data description can be found in section 2, the limitations in section 3, the findings in section 4 and the conclusion is in section 5.

1.1 Research questions

The analysis has been divided into three parts; climactic condition, bushfires and the relationship between climatic condition and bushfires.

Climactic condition:

  • Where does rainfall occur most in Australia? (NEW Q)
  • What is the hottest climate in Australia? (NEW Q)
  • How has global warming impacted temperatures in Australia? (NEW Q)

Bushfires:

  • In which months are bushfires burning?
  • Where are the bushfires burning?

Relationship between climactic condition and bushfires:

  • How does temperature and rainfall affect the number of bushfires?
  • What is the correlation between climactic conditions and bushfires? (NEW Q)

2 Data description

This section mainly introduces the data, data sources and data description.
There are three data sets used on this analysis, and the cleaned data is obtained from GitHub tidytuesday.

2.1 Australian Fire Data

The Australian fire data has been extracted from the MODIS fire product collection at NASA (NASA 2020). The fire data is collected every five minutes and there are 5101817 observations from 2000-11-01 to 2020-01-05. All the variables in the dataset have been presented in the table below. The variables predominantly used in this analysis are; latitude, longitude and acquisition date.

Table 2.1: Australia fire data
Variable Description
latitude Center of 1km fire pixel but not necessarily the actual location of the fire as one or more fires can be detected within the 1km pixel.
longitude Center of 1km fire pixel but not necessarily the actual location of the fire as one or more fires can be detected within the 1km pixel.
brightness Channel 21/22 brightness temperature of the fire pixel measured in Kelvin.
scan The algorithm produces 1km fire pixels but MODIS pixels get bigger toward the edge of scan. Scan and track reflect actual pixel size.
track The algorithm produces 1km fire pixels but MODIS pixels get bigger toward the edge of scan. Scan and track reflect actual pixel size.
acq_date Date of MODIS acquisition.
act_time Time of acquisition/overpass of the satellite (in UTC).
satellite A = Aqua and T = Terra.
confidence This value is based on a collection of intermediate algorithm quantities used in the detection process. It is intended to help users gauge the quality of individual hotspot/fire pixels. Confidence estimates range between 0 and 100% and are assigned one of the three fire classes (low-confidence fire, nominal-confidence fire, or high-confidence fire).
version Version identifies the collection (e.g. MODIS Collection 6) and source of data processing: Near Real-Time (NRT suffix added to collection) or Standard Processing (collection only). ‘6.0NRT’ - Collection 6 NRT processing.’6.0’ - Collection 6 Standard processing. Find out more on collections and on the differences between FIRMS data sourced from LANCE FIRMS and University of Maryland.
dbright_t31 Channel 31 brightness temperature of the fire pixel measured in Kelvin.
frq Depicts the pixel-integrated fire radiative power in MW (megawatts).
day_night D = Daytime, N = Nighttime

2.2 Climate data

The climate data was extracted from the Australian Bureau of Meterology (BoM). The Bureau of Meterology is the weather station that measures rainfall, wind, temperature, etc.

The cleaned Rainfall data was obtained from from GitHub tidytuesday. It has rainfall for six Australian cities, namely; Perth, Adelaide, Melbourne, Sydney, Brisbane and Canberra. It contains more than 230,000 observations and has been collected from 1858-01-01 to 2020-01-06. There was a few missing values for Brisbane and Canberra. To maintain the integrity of the data the missing values have been added from the source website and cleaned so that the data is complete for Canberra (for 1968-01-01 to 2017-12-31) and Brisbane (for 1893-01-01 to 1998-12-31). Therefore, the rainfall dataset used in this analysis is a combination of the two above.

The temperature data has been retrieved from two sources. The first source for temperature data was GitHub tidytuesday, and it has been collected from 1910-01-01 to 2019-05-31. Since the dates for the cleaned temperature data was not consistent with the fire data, a second source was used. The second source for the temperature data was source website. This data was cleaned to obtain temperatures from 2019-06-01 to 2020-01-05. Both datasets were merged to produce the final dataset used for the analysis. The final dataset has around 530,000 observations taken from 1910-01-01 to 2020-01-05. The seven weather stations chosen were based on the seven Australian cities; Perth, Adelaide, Melbourne, Sydney, Brisbane, Port Lincoln and Canberra.

The structure of the climate data is presented in the table below. The year, city name and rainfall variables were mainly used for from the rainfall dataset and the date, temperature and temperature type were predominantly used from the temperature dataset.

Table 2.2: Temperature data
Variable Class Description
city_name character City Name
date double Date
temperature double Temperature in Celsius
temp_type character Temperature type (min/max daily)
site_name character Actual site/weather station
Table 2.3: Rainfall data
Variable Class Description
station_code character Station Code
city_name character City Name
year double year
month character month
day character day
rainfall double Trainfall in millimeters
period double how many days was it collected across
quality character Certified quality or not
lat double latitude
long double longitude
station_name character Station Name

3 Limitations of analysis

The main limitations of the dataset are concerned with the fire and climate data.

There is no regional division in the fire data. Therefore, the data has been assigned to a state or reigion based on its longitude and latitude, which may lead to location bias in the analysis. (borders for NSW an VIc not easily distinguished)

The rainfall and temperature data has been recorded for some major cities, and the sample is relatively small when compared to the fire data. Therefore, the correlation between the fire and climate data in section 4 may not be precisely accurate as the temperature could be for a city which may be hundreds of kilometres away from the fire it has been correlated to. This may cause a deviation in the results.

4 Analysis and findings

4.1 Climate Conditions

4.1.1 Where does rainfall occur the most in Australia?

Average monthly rainfall for 2019

Figure 4.1: Average monthly rainfall for 2019

Figure 4.2: Average monthly rainfall from 2001

Figure 4.1 shows the average monthly rainfall in 2019 for Adelaide, Brisbane, Canberra, Melbourne, Perth and Sydney. From this chart it is evident that Sydney had the most rain throuhgout 2019. The amount of rainfall in each state varies from month to month, which is likely due to each state having slightly different weather patterns and seasons.

4.1.2 Where is the hottest climate in Australia?

## `summarise()` regrouping output by 'year' (override with `.groups` argument)
Daily Average Temperature for each State Overtime

Figure 4.3: Daily Average Temperature for each State Overtime

Figure 4.3 is a heatmap displaying the daily average temperature for each cities. It is clear that the warmest cities are Brisbane and Perth on average. Alarmingly, each area/city is increasing in temperature overtime, which is clear evidence of global warming.

4.1.3 How has climate change impacted weather patterns in Australia?

The plot for the difference between the average temperature of 1961-1990 (as baseline) and the annual average temperature for each year from 1910 to 2019, calculated by daily maximum temperature

Figure 4.4: The plot for the difference between the average temperature of 1961-1990 (as baseline) and the annual average temperature for each year from 1910 to 2019, calculated by daily maximum temperature

Figure 4.4 shows the average temperature difference in Australia from 1910 to 2019. It is compared against the average temperature from 1961-1990. There is a clear upward trend, thus indicating that the average temperature in Australia is increasing. The figure highlights the average temperature for 2019, as it was the highest average temperature on record.

Annual rainfall difference 1961-1990 (as baseline) and the annual average rainfall for each year from 1910 to 2019, calculated by daily maximum rainfall

Figure 4.5: Annual rainfall difference 1961-1990 (as baseline) and the annual average rainfall for each year from 1910 to 2019, calculated by daily maximum rainfall

Figure 4.5 shows the average rainfall difference in Australia from 1910 to 2019. The calculated average rainfall for a given year is compared against the average rainfall from 1961-1990. From 1961-1990 the average rainfall is variable, which is likely a result of the natural weather patterns. However, after 1990 there is a clear downward trend, indicating that each year there is less rainfall than the 1961-1990 average. Interestingly enough, 2019 had one of biggest differences and thus had considerably less rainfall than other years.

Both figures indicate that climate change has lead to Australia becoming a hotter and drier continent. As mentioned above, this is the ‘perfect storm’ for bushfires.

4.2 Bushfires

4.2.1 When are the fires burning?

Figure 4.6: Yearly Australian fires from 2001 until 2020

Figure 4.6 is a line plot showing the number of bushfires for each year, ordered according to month. The trend line for 2019 is clearly different to other years, which is due to the bushfires in Victoria and New South Wales. These were the biggest bushfires since the European Settlement (Nolan et al. 2020), hence the clear upward trend from September 2019.

In addition, there were more bushfires in September 2011 than in any other year in the same period. According to Blanchi et al. (2014), these bushfires were the result of a combination of low rainfall and strong winds.

The overall trend line shows that bushfires are more likely to occur from August to November. This is the winter-spring period in Australia. These months are typically drier in most parts of Australia as seen in Figure 4.1, which can lead to the ignition of forest fuels (Sullivan et al. 2012) and thus bushfires.

4.2.2 Where are the bushfires burning?

Figure 4.7: The number of bushfire in different states or regions in the past 20 years

Most bushfires occur in the Northern Territory and Queensland, as seen in Figure 4.7. This is likely to both areas being prone to drought and have high concentrations of vegetation.

Figure 4.7 shows that the number of fires in Victoria and New South Wales is relatively stable, with the exception of 2019. As mentioned above, this is due to the catastrophic bushfires that took place over this period. Additionally, figure 4.8 shows the distribution of bushfires from the end of 2019 to the beginning of 2020. It predominantly captures Victoria and NSW, due to the extensive fires in 2019.

Figure 4.8: The distribution point for Australian bushfires from 2019-12-29 to 2020-01-05(the darker the color, the more serious the fire)

4.3 How does temperature and rainfall affect the number of bushfires?

4.3.1 Positive association with temperature and bushfires

## `summarise()` regrouping output by 'Year' (override with `.groups` argument)
## `summarise()` regrouping output by 'year' (override with `.groups` argument)
## `summarise()` ungrouping output (override with `.groups` argument)

Figure 4.9: The Trend of Temperature and Bushfires

Figure 4.10: The Average Temperature Difference and Average Bushfires during 2000 until 2020

Figure ?? shows the average temperature in 2019 is significantly higher than other years. A strong positive Indian Ocean Dipole (IOD) phenomenon was the culprit(Harris and Lucas 2019), and contributed to very high temperatures and low rainfall across Australia(Meteorology 2019), beginning in May 2019 and lasting until the end of the year.

Indian Ocean Dipole: positive phase. Source from Australian BOM.

Figure 4.11: Indian Ocean Dipole: positive phase. Source from Australian BOM.

How big is the difference between the temperature in 2019 and other years?

Figure 4.4 shows the average annual temperature from 1961 to 1990. From 1910 to mid-1950, the annual average temperature for almost all years was lower than the baseline. Since the mid-1950s, the average annual temperature has been higher than the baseline. Although there has been the odd year, these are outliers.

Moreover, it is worth noting that 2019 is significantly higher than in other years, with a difference of above 1.5°C. As Meteorology (2019) states, 2019 was Australia’s warmest year on record, surpassing the previous record of +1.35°C in 2013.

Evidently, the rising annual temperature has been a warning that Global warming is a growing problem. Ultimately, global warming leads to many natural disasters, bushfires included.

Figure ?? shows the total number of fires in Australia from 2001 to 2019. Compared with Figure ??, the total number of fires is fairly consistent??. When the temperature is high, there are generally more fires.

The most fires occurred in 2012 with 474,964 fires, however, the average temperature was relatively low, at 22.57°C. These bushfires can be explained by a lightening strike(Dowdy and Mills 2012).

Ultimately, the analysis above demonstrates that a high temperature is closely related to bushfires. When the temperature is high and climate is dry, forest fuels are more likely to catch alight(G. J. van Oldenborgh et al. 2020).

4.3.2 Negative association with rain and bushfires

## NULL

Figure 4.12: The Trend of Rainfall and Bushfires

Figure 4.13: Annual Rainfall and Bushfires

Australia’s location means that rainfall is highly variable. It is strongly influenced by the global climate system phenomena such as El Niño, La Niña, and IOD.

Figure ?? shows how the annual rainfall has changed over time in Australia. In 1994 there was a severe drought, influenced by the El Niño weather pattern. This was the fifth year of drought for some parts of Australia(Nicholls 2004). In 2010, the annual rainfall was the highest in 20 years. This was a result of a La Niña event.

2019 was the driest year in the last 20 years. One reason is that the strongest positive IOD reached the highest values on record across 60 years, and due to the frequency and influence of sea-surface temperature changes, the El Nino-Southern Oscillation(ENSO) is neutral throughout the year??????

As a result of the positive IOD, the Walker Circulation was severely curtailed and there was abnormal easterly winds in the Indian Ocean. These winds meant cloud cover in Australia was swept away, dramatically reducing rainfall across Australia(Hughes 2003).

Figure ?? shows there has been a consistent average number of fires (250,000) in the last 20 years. The number of fires in 2010 was the lowest, which may be due to the high rainfall as a result of the La Niña event.

From Figure 4.5, the annual average rainfall changed cross-time compare with the average baseline. Over the last 20 years, the rainfall has consistently been below average, with the occasional outlier in 2010 and 2011.

It is evident that the low rainfall in 2019 was a contributor to the high number of bushfires.

4.3.3 The correlation between climatic conditions and bushfires

The correlation between rainfall, temperature and bushfires

(#fig:correlation_color)The correlation between rainfall, temperature and bushfires

From Figure @ref(fig:correlation_color), it seems that the rainfall and bushfires have small negative correlation, which is only 30%. This means that if the average annual precipitation increases, the probability of bushfires may decrease with 30% of the increase in precipitation.

As for temperature and bushfires, the correlation between these two variables is only 0.13 which is too close to 0. Therefore, Our group believes that there is no correlation between temperature and occurrence of bushfires in the past 20 years. More data may be needed later to draw the different conclusion.

Finally, let’s observe the precipitation and temperature. The correlation between them is -0.47, which is the largest absolute value among the three correlations, which indicates that their correlation degree is the highest. The correlation coefficient of 0.47 means that if the temperature rises, the precipitation will decrease to 47% of the range that rise. In my opinion, this is because temperatures will impact the rate of evaluation, with higher temperatures leading to drier weather and faster air moisture loss. (Hausfather 2018).

In conclusion, from corrplot it is clear that the main climatic conditions to affect bushfires is the decreasing of the rainfall. What surprided our group was that the direct effect of temperature rise on Australian bushfires occurrence was almost 0. However, due to the temperature will affect the precipitation and thus the occurrence of bushfires, so people should pay attention to the greenhouse effect and other issues. The upper graphs show that 2019 is the hot and drive year in Australia, which also has the most massive bushfires occurred. That discovery confirms the accuracy of our corrplot again. The increase of temperature leads to the decrease of rainfall have been a primary driver of increases in bushfires area of Australia. (G. J. van Oldenborgh et al. 2020)

5 Conclusions

This analysis has looked at when, where and what climactic condition bushfires are likely to thrive in. Bushfires are common from August to November, and are frequently found in the Northern Territory and Western Australia. High temperature and drought are the two biggest contributors to bushfires(G. J. van Oldenborgh et al. 2020). This is consistent with data from 2019, which discovered that it was the hottest and driest year on record with the highest number of bushfires.

Acknowlegments

The authors would like to thank all the contributors to the following R package: Wickham et al. (2019), Wickham (2016), Wickham, Hester, and Francois (2018), Cheng, Karambelkar, and Xie (2019), Ryan and Ulrich (2020), Müller (2017), R Core Team (2020), Arnold (2019), Wickham et al. (2020), Vanderkam et al. (2018), Grolemund and Wickham (2011), Schloerke et al. (2020), Rudis (2020).

6 References

7 Year section

7.1 VEGA"S PART (ADDITIONAL TO MASTER)

7.1.1 In which months are bushfires burning?

When the trend in other years (i.e. 2011 and 2012) seems to reach a peak in August and October respectively, the trend of bushfire in 2019 is still growing until December 2019. Unfortunately, we do not have additional data to see the continuity of the 2019 trend line in 2020.

According to the figure, the most severe bushfire occurs in 2011 with the number of spots reach 120.348 fire spots. We will look further into 2011 in the figure below.

–>

7.1.2 Where are the bushfires burning?

Furthermore, we also draw a line that shows national average number of bushfires during 20 years period. Hence, we can see which states/territories that have more bushfires compared to national average. In addition, we are interested to see particular year of 2011 (the worst year) and 2019 (the different pattern).

Figure 4.7 shows us that, in general, bushfires are mostly occured in the Northern Territory(NT) and West Australia (WA), possibly because both areas are more likely prone to drought (we will check this assumption in the following section) and comprises of large forests area according to Australia’s Department of Agriculture.

In 2011, NT has very high number of bushfires compared to other states. However, latest massive bushfires in 2019 is mostly occured in area comprises of Victoria and New South Wales. Figure 4.8 shows fire spots from the latest bushfires (December 2019 to early January 2020), and we can find that it is mainly concentrated in Victoria and NSW.

Based on the analysis, we will look at the causes of bushfires in Australia. Are there association between climate conditions to the bushfire? We will explore the story of climate conditions (temperature and rainfall level) in section 4.3.1 and 4.3.2.

7.2 How does temperature and rainfall affect the number of bushfires?

7.2.1 Positive association with temperature and bushfires

This section will be focus on finding clues to answer research question of whether higher temperature contributes to higher frequency of bushfires occurences. Our base assumption is that these two variables have a positive association which means that if the temperature increase the number of fire spots should also increase and vice versa.

In order to see the association between temperature and bushfire, we will produce a plot that compares these two variables side by side. To make a more visible pattern, temperature variable will be plotted as a temperature difference from the average temperature during period of observations (2000 - 2020). Mathematically, the temperature difference is formulated as below:

\[ Average\;Temperature\;Difference_n = Annual\;Temperature_n\; - Average\;Temperature\;in\;Period\;of\;Observations\]

On the other hand, since the average number of bushfires (fire spots) has different scale with the temperature, we will scale down the annual average bushfires using log 10 and square root. By doing so, we can put the two variables in a single plot to see the patterns of each variables.

As we can see from figure 4.4, the annual temperature trend seems to rise since 2000 until 2019. Likewise, the annual bushfires also slighlty increase during observed years. Therefore, we can see a positive association between these two variables, although, the association seems not really strong.

As an example, we can see that in 2011 the average bushfires rise sharply (marked by black dashed line) while the temperature only a bit increase from 2010. Another example of the weak association is also shown in the trends during 2015 until 2020. In that period, average bushfires seems to slighlty increase whereas the temperature actually rises significantly. Further analysis on correlation between temperature and bushfires will be discussed in the last section of this paper.

7.2.2 Negative association with rain and bushfires

Has the less rainfall affected the more bushfires in Australia? Or has the more rainfall influenced the less bushfires? Because of the location of Australia, the rainfall in Australia is highly variable, which is strongly influenced by global climate system phenomena such as El Niño, La Niña, and IOD. Despite this large natural variability, the potential long-term trends are evident in some regions, even effecting the local rainfall.
In this section, we will explore the relationship between rainfall and bushfires, is it positive or negative?

As a first step, we observe annual rainfall level that retrieved from the rainfall dataset. The interval of observation is ranging from 2000 until 2020.

Similarly with our previous job when comparing the trend of average temperature to bushfires, we will rescale down the annual average fire spots using log 10 in order to put its trend side by side with average rainfall trend in the same plot.

As shown in figure 4.13, we can see that, in general, annual rainfall has a decreasing trend during the observed period (2000 - 2020), while bushfires has a slightly rising trend. Hence, we can assume that rainfall and bushfires are negatively associated. This contrast pattern is quite obvious particularly in 2011 where average fire spots soar as the average rainfall level drops from 4.9 mm in 2010 to 4.15 mm in 2011. After 2015, bushfires is slightly increasing until 2019 as rainfall level continues to drop to its lowest level (2.16 mm) in the same year.

To complete this analysis, we will check the exact coefficient number of correlation between variable rainfall and bushfires in the following section.

Arnold, Jeffrey B. 2019. Ggthemes: Extra Themes, Scales and Geoms for ’Ggplot2’. https://CRAN.R-project.org/package=ggthemes.

Australia, Geoscience. 2020. “Bushfire.” https://www.ga.gov.au/scientific-topics/community-safety/bushfire.

Blanchi, Raphaele, Justin Leonard, Katharine Haynes, Kimberley Opie, Melissa James, and Felipe Dimer de Oliveira. 2014. “Environmental Circumstances Surrounding Bushfire Fatalities in Australia 1901–2011.” Environmental Science & Policy 37: 192–203.

Cheng, Joe, Bhaskar Karambelkar, and Yihui Xie. 2019. Leaflet: Create Interactive Web Maps with the Javascript ’Leaflet’ Library. https://CRAN.R-project.org/package=leaflet.

Dowdy, Andrew J, and Graham A Mills. 2012. “Characteristics of Lightning-Attributed Wildland Fires in South-East Australia.” International Journal of Wildland Fire 21 (5): 521–24.

Grolemund, Garrett, and Hadley Wickham. 2011. “Dates and Times Made Easy with lubridate.” Journal of Statistical Software 40 (3): 1–25. http://www.jstatsoft.org/v40/i03/.

Harris, Sarah, and Chris Lucas. 2019. “Understanding the Variability of Australian Fire Weather Between 1973 and 2017.” PloS One 14 (9): e0222328.

Hausfather, Zeke. 2018. “Explainer: What Climate Models Tell Us About Future Rainfall.” Carbon Brief.

Hughes, Lesley. 2003. “Climate Change and Australia: Trends, Projections and Impacts.” Austral Ecology 28 (4): 423–43.

Meteorology, Bureau of. 2019. “Annual Australian Climate Statement 2019.”

Müller, Kirill. 2017. Here: A Simpler Way to Find Your Files. https://CRAN.R-project.org/package=here.

NASA. 2020. “NASA Fires of Australia.”

Nicholls, Neville. 2004. “The Changing Nature of Australian Droughts.” Climatic Change 63 (3): 323–36.

Nolan, Rachael H, Matthias M Boer, Luke Collins, Vı́ctor Resco de Dios, Hamish Clarke, Meaghan Jenkins, Belinda Kenny, and Ross A Bradstock. 2020. “Causes and Consequences of Eastern Australia’s 2019–20 Season of Mega-Fires.” Global Change Biology 26 (3): 1039–41.

Oldenborgh, Geert Jan van, Folmer Krikken, Sophie Lewis, Nicholas J Leach, Flavio Lehner, Kate R Saunders, Michiel van Weele, et al. 2020. “Attribution of the Australian Bushfire Risk to Anthropogenic Climate Change.” Natural Hazards and Earth System Sciences Discussions, 1–46.

Oldenborgh, Geert Jan van, Folmer Krikken, Sophie Lewis, Nicholas J Leach, Flavio Lehner, Kate R Saunders, Michiel van Weele, et al. 2020. “Attribution of the Australian Bushfire Risk to Anthropogenic Climate Change.” Natural Hazards and Earth System Sciences Discussions, 1–46.

R Core Team. 2020. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing. https://www.R-project.org/.

Rudis, Bob. 2020. Hrbrthemes: Additional Themes, Theme Components and Utilities for ’Ggplot2’. https://CRAN.R-project.org/package=hrbrthemes.

Ryan, Jeffrey A., and Joshua M. Ulrich. 2020. Xts: EXtensible Time Series. https://CRAN.R-project.org/package=xts.

Schloerke, Barret, Di Cook, Joseph Larmarange, Francois Briatte, Moritz Marbach, Edwin Thoen, Amos Elberg, and Jason Crowley. 2020. GGally: Extension to ’Ggplot2’. https://CRAN.R-project.org/package=GGally.

Sullivan, Andrew L, W Lachie McCaw, Miguel G Cruz, Stuart Matthews, and Peter F Ellis. 2012. “Fuel, Fire Weather and Fire Behaviour in Australian Ecosystems.” Flammable Australia: Fire Regimes, Biodiversity and Ecosystems in a Changing World, 51–77.

Vanderkam, Dan, JJ Allaire, Jonathan Owen, Daniel Gromer, and Benoit Thieurmel. 2018. Dygraphs: Interface to ’Dygraphs’ Interactive Time Series Charting Library. https://CRAN.R-project.org/package=dygraphs.

Wickham, Hadley. 2016. Ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York. https://ggplot2.tidyverse.org.

Wickham, Hadley, Mara Averick, Jennifer Bryan, Winston Chang, Lucy D’Agostino McGowan, Romain François, Garrett Grolemund, et al. 2019. “Welcome to the tidyverse.” Journal of Open Source Software 4 (43): 1686. https://doi.org/10.21105/joss.01686.

Wickham, Hadley, Romain François, Lionel Henry, and Kirill Müller. 2020. Dplyr: A Grammar of Data Manipulation. https://CRAN.R-project.org/package=dplyr.

Wickham, Hadley, Jim Hester, and Romain Francois. 2018. Readr: Read Rectangular Text Data. https://CRAN.R-project.org/package=readr.